Molecular Nanostructures on the Surface of a cL^^-Superconductor 
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We study molecular nanostructures on the surface of a d x 2_ y 2-w&ve superconductor. We show 
that the interplay between the molecular nanostructure's internal excitation spectrum and quantum 
interference of the scattered host electrons leads to a series of novel effects in the local density of 
states. We demonstrate that these effects give insight into both the nature of the superconducting 
pairing correlations and of the intermolecular interactions. 
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Nanostructures provide the intriguing possibility of 
manipulating in a controlled way the electronic struc- 
ture of the host system they reside on This level of 
control opens new venues for studying the complex elec- 
tronic structure of many strongly correlated electron sys- 
tems. For example, it was suggested that nanostructures 
can elucidate the nature of superconducting correlations 
in conventional Q and unconventional superconductors 
0. Nanostructures formed from more complex build- 
ing blocks, such as molecules with internal degrees of 
freedom, provide a novel way of manipulating and simul- 
taneously gaining insight into the electronic structure of 
complex systems. The internal vibrational and rotational 
excitations of single molecules have been intensively stud- 
ied by inelastic tunneling spectroscopy 0, 0| . Their in- 
teraction with and effect on the electronic structure of 
a metallic surface was recently investigated by Gross et 
al. 0. 

In this Letter we study nanostructures composed of 
two molecules that reside on the surface of an unconven- 
tional d x 2_ y 2 -superconductor that represents the family 
of high-temperature superconductors (HTSC). We show 
that the interplay between the molecular nanostructure's 
internal excitation spectrum and quantum interference 
of the scattered host electrons leads to a series of novel 
effects in the local density of states (LDOS) of the su- 
perconductor. Our main results are threefold. First, we 
identify three types of intermolecular interactions and 
show that each of them leads to qualitatively distinguish- 
able features in the LDOS. This, in turn, allows one 
to identify the nature of intermolecular interaction from 
experimental measurements of the LDOS. Second, we 
demonstrate that the LDOS changes with the molecules' 
distance and orientation relative to the underlying host 
lattice. This effect permits us to probe the spatial depen- 
dence of superconducting correlations. Third, we show 
that by exciting specific energy levels of the molecular 
structure, the LDOS can be manipulated in a controlled 
manner. These results intricately relate the study of 
strongly correlated electron systems with the further de- 
velopment of molecular electronics. 

A nanostructure consisting of N molecules, repre- 
sented by local bosonic modes Jjj , on the surface of a two- 



dimensional d x 2_ y 2 superconductor possesses the Hamil- 
tonian H = H e + + Hi nt where 

H e = ^ e k4,<r C k^ + AkC M c -k,l + h - c - » 
k, CT =T4 k 

H b =u> o y^ y a\. a Fi + J ^ [c4 f , Q r< , aj, , a r ,.] , 



(1) 



H e and Hb describe the unperturbed superconductor 
and nanostructure, respectively, and Hi nt represents 
the interaction between them. c\a^ are the fermionic 
and bosonic creation operators, respectively. A k = 
A (cos k x cos hy) /2 is the <i a ,2_y2-wave superconducting 
gap and e k = — 2t (cos k x + cos k y ) — At' cos k x cos k y — fi 
is the normal state tight-binding dispersion, ujq is the 
characteristic frequency of the N modes located at i"; 
(i = 1, .., N). ^ is a quadratic functional in the bosonic 
operators and describes the interaction between modes 
at sites and Tj with strength J (three specific forms of 
*Sf are discussed below), g is the boson- fermion scatter- 
ing vertex. The fermionic Greens function of the unper- 
turbed (clean) system in Nambu notation is 



G 1 (k, iu n ) = [iu n T - e k r 3 ] <r + A k T 2 cr 2 



(2) 



where <Xj and 7j are the Pauli matrices in spin and Nambu 
space, respectively. Diagonalizing H b via a Bogoliubov 
transformation to new operators bj,bi (1 = 1,..., N), one 
obtains H b = Slib\bi and 



H\, 



where g£l is a site and mode dependent interaction ver- 
tex, and ill is the bosonic energy spectrum. The full 
fermionic Green's function is now given by 

G(r, r>„) = G (r, r',oj n ) + G (r, r i; w„) 

xS(r l ,r J , w„)G (rj,r', uj n ) , (4) 



i.l.a 



(b\ + &,) 4 



(3) 
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where E is the full fermionic self-energy obtained from 
the Bethe-Salpeter equation 

E^r,-,^) = £o(rj, rj,u) n ) + Epfa, r p , u n ) 

xG (r p ,r (? ,w rl )S(r g ,r i ,w Il ) (5) 
with p,q = 1, N and 

£ (ri,rj,u n ) = T r 3 ffff G (r» , r j ,u) n + v m ) 

m,l 

xD^rsgi 1 ] . (6) 

The retarded form of the bosonic propagator is given by 

Df-(w) = [uj + iY(uo) - VLiY 1 - [uj + W(u) + Qir 1 . (7) 

The lowest order vertex corrections scale as Sg/g ~ 
(gkp/4:Vp) 2 F where F is a function of OCV) for bosonic 
and fermionic frequencies smaller than A |8j . For the pa- 
rameter range considered below, gkp/AvF <C 1 and ver- 
tex corrections can thus be neglected. The interaction, 
Hi n t, not only leads to changes in the fermionic LDOS, 
but also affects the bosonic excitation spectrum twofold. 
First, Hi n t shifts the unperturbed bosonic frequencies, a 
shift that we take to be included in the effective values for 
considered below. Second, the bosonic modes acquire 
a finite lifetime, V(lo) ^ 0. The lowest order bosonic self- 
energy correction yields T(ui) = g 2 u 3 / (6nvpv\), where 
vf («a) is the Fermi (superconducting) velocity at the 
nodal points perpendicular (parallel) to the Fermi sur- 
face. Since for all cases considered below T(lu ) -C Wo, 
which requires g <C 50Ao, T(o;) can be neglected and 
we take for simplicity T(oj) — + . Finally, the LDOS, 
N = An + M2 with A u (y,uj) = -21mGu(r, uj + iS) is 
obtained numerically from Eqs.Q with 5 = 0.2 meV. 

We showed in Ref. that a single local mode on the 
surface of a d x i _ y 2 -superconductor is pair-breaking and 
induces a fermionic resonance state, whose spectroscopic 
signature are two peaks in the LDOS. Introducing a sec- 
ond mode leads to two new physical effects that introduce 
qualitative changes in the LDOS. First, due to quantum 
interference of the scattered electrons, the calculation of 
the LDOS involves not only the local Greens function, 
Go(rj,r,-), and self-energy, £(ri,rj), with I - .; = rj, but 
also the non-local ones with ^ Yj. As a result, the 
LDOS is determined by the distance and orientation of 
the modes relative to the lattice, which is described by 
Ar = Yi — Yj — (n x ,n y ) (we set the lattice constant 
ao = !)• Second, the intermolecular interaction ^ leads 
to characteristic features in the LDOS that, as we show 
below, depend on the specific form of ^. 

In order to distinguish between the effects of quan- 
tum interference, superconducting correlations and inter- 
molecular interaction on the LDOS, we consider first a 
simplified electronic bandstructure with t = 300 meV, 



/i = t' = and set J = 0. In this case, superconduct- 
ing correlations, as represented by the non-vanishing of 
the non-local anomalous Greens function, F^r, and self- 
energy, <&Ar (the off-diagonal elements of Go and S , 
respectively), are present for (n x + rij,)mod2 = 1 (case 
T) and vanish otherwise (case IT). We find that the pres- 
ence of superconducting correlations leads to qualitative 
differences in the LDOS between case I and II that are 
universal, i.e., independent of the specific realization of 
case I or II considered. For definiteness, we consider be- 
low two modes located at ri = (0, 0) and r 2 = (1, 0) (case 
I) and ri = (0, 0) and r 2 = (2, 0) (case II) and present 
the resulting LDOS at T = in Fig. Ufa) and (b). The 




Frequency [Ao] Frequency [Ao] 

FIG. 1: LDOS at R as a function of frequency for loo ~ O.6A0 
and three intermolecular interactions: (a),(b) HO coupling, 
(c),(d) FM coupling, and (e),(f) AFM coupling. Scaling of 
the curves is indicated. 

coupling of the resonance states due to quantum inter- 
ference only (see curves for J = 0) leads to the formation 
of bonding and antibonding states and an energy split- 
ting between them. Only one of these states is located at 
energies smaller than Ao, resulting in a particle-like and 
hole-like peak in the LDOS. Quantum interference leads 
to oscillations in the frequency of the resonance peaks 
when Ar is changed (cf. Figs, ^a) and (b)), similar to 
the case of non-magnetic impurities Next, we discuss 
the effects of three different intermolecular interactions 
on the frequency and shape of the resonance peaks. 

Harmonic Oscillator (HO) Coupling: A coupling of 
the bosonic modes via a quadratic potential is the 
quantum mechanical analog of two harmonic oscilla- 
tors coupled by a spring. In this case ^ = 
(6^ + b ri — bl 2 — & r2 ) with J > 0. After diagonaliza- 
tion, one has fii = ojq, fl2 = ^oV^ + 4A with A = J/loq, 
gill = 9/V2 and gi% = ±g/[V2(l + 4A) 1 /*]. The re- 
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suiting LDOS for several values of J is shown in Fig.QJa) 
and (b) for cases I and II, respectively. With increasing 

(1) 

J, il 2 increases while \gr 13 1 decreases, the electronic scat- 
tering hence becomes weaker and the resonance peaks are 
shifted to higher energies. In case I, the peaks move close 
to ±Ao as J increases from J = to J = oo, an effect 
that is independent of g and directly related to a non-zero 
i^Ar and &Ar- In contrast, in case II, the frequency shift 
is much smaller. For J = oo, only the bosonic in-phase 
mode couples to the fermionic system. This coupling is 
much weaker in case I than in case II due to the super- 
conducting correlations. 

Ferromagnetic (FM) Coupling: A second type of in- 
termolecular interaction arises if each molecule possesses 
(2L + 1) states that are represented by pseudo-spins Li 2, 
whose interaction is given by the anisotropic Heisenberg 
Hamiltonian 

H S pi n — J z LfL 2 + J± (LJZf ~t~ ^1-^2) • (^) 

For J z < and \J Z \ > |J±|, the pseudo-spins are ferro- 
magnetically aligned (i.e., the molecules prefer to be in 
the same state). Performing a Holstein-Primakoff trans- 
formation followed by a large- L expansion up to order 
O(L), we obtain H b in Eq.JU with ujq = —J Z S > 0, 
J = J±S > and * = bljrz + 6j 2 6 ri . After diago- 

nalization, one has fl^a = u>o T J, fc, = Tg/V% and 

5ri^ 2 = g/V%- The operators 6j. , 6 ri represent the exci- 
tations of the pseudo-spin dimcr which scatter fermions 
via Hi n t. For J z < J±, i.e., J > u>o, the system be- 
comes unstable since the ferromagnetic alignment of the 
pseudo-spins is destroyed. 

For small values of g, the resonance peaks move to- 
wards lower energies in cases I and II with increasing 
J (not shown). However, for intermediate g w 30A , 
this behavior changes qualitatively for case I; the res- 
onance frequencies exhibit a non-monotonic dependence 
on J [Fig.njc)], in contrast to case II [Fig.^d)]. For even 
larger g, the resonance peaks shift to higher energies in 
case I, and to lower energies in case II with increasing J. 
A detailed analysis shows that this qualitative difference 
arises from a non-zero F& r and <I>Ar m case I. Note, in 
case II (Fig. ^d)), the width of the resonance peaks in- 
creases when it shifts to lower energies, in contrast to the 
effect expected from a decreasing residual DOS 0] . Here, 
however, the width of the resonance peaks is also deter- 
mined by the imaginary part of the self-energy which van- 
ishes for |Q| < Since Sli decreases with increasing J, 
the resonance state becomes more strongly damped once 
ill < w res . In Fig. njd), we have ill — ujq > uj res for 
J = 0, but rii < uj res for J = 2ujq/3 (Qi = 0.2Ao) and 
J = luq (ill = 0), resulting in an increased peaks' width 
for the latter two values of J. 

Antiferromagnetic (AFM) Coupling: For antiferro- 
magnetic alignment of the pseudo-spins (J 2 > 0, \J Z \ > 
J± > in Eq.JSJl), the molecules prefer to be in "op- 



posite" states. Performing a Holstein-Primakoff trans- 
formation followed by a large-L expansion, one obtains 
* = ^ri^r, + b T1 b T2 resulting in two degenerate energy 
dispersions 0,i j2 = c^oVl — A 2 with A = J/u>o, gll' 2 ^ — 
±gu-L, 2 , gll^ = Tgu 2 ,i and u\ 2 = [1/Vl - A 2 ± l]/2. 
For J > cJo ; i.e., \J Z \ < J±, the system becomes unstable 
since the antiferromagnetic alignment of the molecular 
pseudo-spins is destroyed. 

(1 2) 

With increasing J, S2 1>2 decreases while g^ 2 ' increases, 

thus leading to an increase in the scattering strength. 

As a result, the resonance peaks shift monotonically to 

lower energies in case II [see Fig. [Iff)], while their width 

decreases since uj res < O12 for all J. In the limit J — > luq, 
(12). 

g r2 diverge, and the scattering becomes unitary. Since 
the local and non-local Go and E vanish for to = 0, the 
resulting resonance peaks are located at zero energy and 
the LDOS vanishes at the molecules' sites. In contrast, in 
case I [Fig. me)], the resonance peaks first shift to lower 
energies with increasing J, but eventually move back to 
higher energies and broaden. This behavior arises from 
the nonvanishing of the non-local Greens function and 
self-energy for w = 0. 

At T = 0, the imaginary part of the normal self- 
energy (the diagonal element of So) possesses logarith- 
mic divergences that lead to dips in the LDOS at ±u>f 2 
(ujf 2 = A Q ± fh j2 ) ISE3- At T ^ 0, the positive en- 
ergy levels of the bosonic modes, represented by the first 
term in Eq.J7J), become populated, opening a new chan- 
nel for fermionic scattering. As a result, Eo acquires 
logarithmic divergences at ±u>i 2 - This leads to either 
peaks or dips in the LDOS, depending on whether lj~ 
is smaller or larger than the frequency of the resonance 
peaks, ±u res . The LDOS for T ^ is shown in Fig. H 
In Figs.^fa) and (b), we plot the LDOS for HO coupling 
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FIG. 2: Temperature dependence of the LDOS for u>o = 
0.6A and J = 2lj /3: (a), (b) HO-coupling, (c) FM-coupling, 
and (d) AFM-coupling. 

in case II for g = I5Ao and g = 24Ao, respectively. For 
g = 15Aq, lui 2 < Lo res and the logarithmic divergences 
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FIG. 3: LDOS at T = for lu = 0.6A and g = 15 A . 
(a),(b) The positive energy branches of the bosonic modes 
are selectively populated for HO-coupling with n^ 2 = 1 and 
J = uo/3. (c), (d) LDOS for a HTSC band structure and (c) 
HO-coupling, and (d) AFM-coupling. 

lead to four peaks (indicated by arrows) in the LDOS. 
Their intensity is governed by n^f^) and thus larger at 
±u>^ than at ±to 2 ■ Simultaneously, opening of the new 
scattering channels leads to an increase in the imaginary 
part of the self-energy, resulting in a shift of the reso- 
nance peaks to higher energies and an increase in their 
width. In contrast, for g = 24 Ao, loi ~ w res , leading 
to a splitting of the resonance peaks and to a dip in the 
LDOS. The second set of divergences at uj 2 induces two 
peaks in the LDOS, which due to the overall increase in 
the LDOS at lower energies are barely perceptible. The 
LDOS for FM coupling and AFM coupling in case II is 
shown in Fig. [2(c) and (d), respectively. For FM cou- 
pling, the LDOS exhibits two peaks and two dips since 
lu± < Lu res < uj 2 ■ In contrast, for AFM coupling, the 
two modes are degenerate and since w-f 2 > ui res , the log- 
arithmic divergences result in a single set of dips. The 
LDOS in case I (not shown) exhibits qualitatively sim- 
ilar behavior to the ones discussed above for all three 
couplings. Note that the frequencies of the features aris- 
ing from the logarithmic divergences are independent of 
g, which thus permits a direct measurement of ±tof 2 . 

The degeneracy of the modes for AFM coupling results 
in a LDOS with only two features at ±a>j~ and ±o;J~, re- 
spectively, in qualitative contrast to the LDOS for FM 
and HO coupling. The latter two couplings, however, 
can be distinguished by comparing ujq obtained from 
the LDOS near a single molecule J8J with tof 2 . Thus 
the nature of the intermolecular interaction, similar to 
a molecule's internal structure 0, can be identified ex- 
perimentally from the specific features it induces in the 
LDOS. 

Changing the population, nf 2 , of the energy levels 
±01,2, for example by optical means [Til ], opens intrigu- 
ing venues to manipulate the local electronic structure of 
the host material. In Fig.|3(a),(b), we present the LDOS 



for HO coupling with n l 2 — 1 and n[ 2 ^ 0. A non-zero 
n\ [Fig. 12(a)] induces dips in the LDOS at ±w 1 ~, while an 
increase of n 2 [Fig. CJb)] leads to peaks at ±u> 2 , similar 
to the temperature induced changes discussed above. 

Finally, the use of a bandstructure representative of 
the HTSC with t'/t = -0.4 and n/t = -1.18 12] mixes 
the effects of intermolecular interaction, superconducting 
correlations and quantum interference on the LDOS. For 
example, the J-dependence of the LDOS for HO coupling 
[Fig-I3 c )] is similar to the one shown in Fig. da), while 
the LDOS for AFM coupling [Fig. 0(d)] shows a depen- 
dence on J that is qualitatively different from that in 

Fig.rnf). 

In conclusion, we study the effects of molecular nanos- 
tructures on the electronic structure of a d a .2_ y 2-wave 
superconductor. We show that different intermolecular 
interactions lead to qualitatively distinguishable signa- 
tures in the LDOS. By changing the population of the 
bosonic excitations, one can manipulate the host's elec- 
tronic structure and gain further insight into the nature 
of unconventional superconducting correlations. 
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